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Abstract. 

Very high energy (> 100 GeV) gamma-ray emission has been detected for the first time from the composite supernova remnant 
G 0.9+0.1 using the H.E.S.S. instrument. The source is detected with a significance of « 13 a, and a photon flux above 200 GeV 
of (5.7 ± 0.7 IMf + 1.2, vs ) x 10~ 12 cm~ 2 s _1 , making it one of the weakest sources ever detected at TeV energies. The photon 
spectrum is compatible with a power law (dN/dE oc E~ T ) with photon index T = 2.40 ± 0.1 1 IM ± 0.20, VJ . The gamma-ray 
emission appears to originate in the plerionic core of the remnant, rather than the shell, and can be plausibly explained as 
inverse Compton scattering of relativistic electrons. 
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1. Introduction 



2. H.E.S.S. Observations and Results 



G 0.9+0.1 is a well known composite supernova rem- 
nant, recognised as such from its radio morphology 
(Helfand & Becker fl987t . It exhibits a bright compact core 
(~2' across) surrounded by an 8' diameter shell. The radio 
spectrum of the core is significantly harder (a « 0. 12) than that 
of the shell (a «-0.77) (LaRosa et al. l20"U0l X-ray observa- 
tions of the nebula with BeppoS AX (Mereghetti et al. [1998 1, 
Chandra (Gaensler et al. 12001b and XMM-Newton (Porquet 
et al. 2003 1 have unambiguously identified the core region 
as a pulsar wind nebula (PWN). A plausible candidate for 
the central pulsar is the hard spectrum point source CXOU 
J174722. 8-280915; however no pulsed emission has been de- 
tected. Observations with XMM have revealed a softening of 
the X-ray spectrum with increasing distance from the centre 
of the PWN, a signature of energy loss of electrons within the 
nebula. The location of G 0.9+0.1 in the Galactic Centre (GC) 
region suggests a distance of ~8.5 kpc (Mezger et al. [T996l. 
At this distance the size of the shell implies an age of a few 
thousand years for the SNR. At gamma-ray energies G 0.9+0. 1 
has not been detected previously. It is not associated with any 
EGRET source and the only published observation in the TeV 
domain is an unconstraining upper limit from the HEGRA 
collaboration (flux <4.6 x 10~ 12 cnrV above 4.2 TeV, 
Aharonian et al. 12002b . 

The High Energy Stereoscopic System (H.E.S.S.) is a 
newly completed instrument designed to study astrophysical 
gamma radiation in the energy range 100 GeV - 10 TeV. 
H.E.S.S. consists of four imaging atmospheric Cherenkov tele- 
scopes (Hinton l2004l Bernlohr et al. 120011 Vincent et al. 120031 
Aharonian et al. 2004a i; shower images are required in at least 
two of them to trigger the detector (Funk et al. 2004). Using the 
stereoscopic technique the instrument reaches an angular reso- 
lution of ~0.1° and a point source sensitivity of 1% of the flux 
from the Crab Nebula (5 cr in 25 hours). 

Observations of the GC region were made with the par- 
tially complete H.E.S.S. system in 2003. These data revealed 
emission from the Sgr A region, consistent with the position of 
Sgr A* (Aharonian et al. 2004b I. Within the same field of view 
indications for a signal from the supernova remnant G 0.9+0.1 
were seen at the 4<x level. A relatively deep H.E.S.S. obser- 
vation of this region with the complete telescope system was 
performed in March-September 2004. The results from this 
dataset relevant to G 0.9+0. 1 are presented here. A detailed dis- 
cussion of the implications of the 2004 dataset for the region 
around Sgr A* will be presented elsewhere. 
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The large field of view of the H.E.S.S. cameras (~5°) pro- 
vides good sensitivity for point sources at an angular distance 
up to ~2° from the pointing direction of the telescope system. 
The bulk of observations presented here were taken in wobble 
mode around the position of SgrA*. In this mode 28-minute 
runs are taken pointing 0.5° away from the nominal source 
in alternating directions. Approximately 20% of the data were 
taken in wobble mode around G 0.9+0.1 itself. The observa- 
tions have an average offset of 0.9° from the SNR. The off- 
axis sensitivity of the system derived from Monte Carlo simu- 
lations has been confirmed via observations of the Crab Nebula 
(Aharonian et al. 2004c I. The total live time of the observations 
after run selection is 50 hours, from a total observation time of 
60 hours. The mean zenith angle of observation was 18°, re- 
sulting in an energy threshold of 170 GeV after standard cuts. 




Fig. 1. Gamma-ray point source significance map for the re- 
gion around the GC. The position of G 0.9+0.1 is marked with 
a triangle. SgrA* is marked with a star. The circles show the in- 
tegration regions contributing to the quoted significance at the 
positions of these two objects. The six telescope pointings are 
shown as crosses. 

Standard event reconstruction was applied to these data: 
tail-cuts image cleaning, Hillas parameterisation and stereo- 
scopic direction reconstruction based on the intersection of im- 
age axes (Aharonian et al. 2004d). Cuts on the scaled width and 
length of images (optimised on gamma-ray simulations and off- 
source data) are made to suppress the hadronic background. An 
additional cut on the minimum size of images (200 photoelec- 
trons) is made to further suppress the background and select a 
subset of events with a superior angular resolution (0.07°rms), 
but with an increased energy threshold of 350 GeV. Events with 
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directions reconstructed within an angle 6 of a trial source lo- 
cation are considered on-source where 6 < 0.1°. A ring of 
radius 0.5° and an area 7 times that of the on-source region is 
used to derive a background estimate. The on-source and back- 
ground region counts, together with a normalisation factor for 
the different acceptance of these regions, are used to derive the 
statistical significance of any excess following the likelihood 
method of Li & Ma l fT983l 

Figure ^ shows a significance sky map for the field of 
view of the H.E.S.S. GC observations derived as described 
above. Two sources of very high energy (VHE) gamma-rays 
are clearly visible: the Galactic Centre (HESS J1745-290) and 
a second source (HESS J1747-281) coincident with the SNR 
G 0.9+0.1. The statistical significance of HESS J1745-290 is 
35 <x in this dataset, compared with 1 1 cr in our 2003 data 
(Aharonian et al. 2004b I. This increase is expected (for a steady 
source) from the greater sensitivity of the full 4-telescope array 
with respect to the 2-telescope configuration used in 2003, and 
from the increased exposure time in 2004. 

The H.E.S.S. 2003 dataset showed evidence for gamma-ray 
emission at the position of G 0.9+0. 1 at the 4 cr level. This evi- 
dence is confirmed by the presence of a 13 cr excess at the same 
position in the independent 2004 dataset. The gamma-ray like 
excess is well fit by the point spread function of the instrument, 
yielding a best fit position of I = 0.872° ± 0.005°, b = 0.076° ± 
0.005°, consistent within statistical errors with the position of 
the PWN in GO.9+0.1 (/ = 0.871°, b = 0.0772°) (Gaensler 
et al. 2001 1. The systematic pointing error is approximately 20" 
in each direction for these data. 

The simulated point-spread function of the instrument 
provides an adequate description of the excess. Assuming 
a radially symmetric Gaussian emission region (p oc 
exp(-# 2 /2cr 2 ource )) a 95% confidence limit on the extension 
0"source < 1 .3' is obtained. For emission from a uniform thin 
shell a limit of 2.2' on the shell radius can be derived at 
the same confidence level. Figure |2] shows a radio map of 
GO.9+0.1 at 90 cm (LaRosa et al. 2000 1 showing the SNR 
shell and the central compact source (colour scale). Contours 
of VHE gamma-ray emission are superimposed along with the 
point spread function of the instrument for comparison. The 
reconstructed gamma-ray source position is marked as a cross. 
Given the close positional coincidence, we identify the VHE 
gamma-ray source HESS J1747-281 with GO.9+0.1. 

Figure shows the reconstructed gamma-ray spectrum of 
GO.9+0.1 using two sets of event selection cuts. A size cut 
of 80 photoelectrons (pe) is applied to extend the spectrum to 
lower energies. This dataset can be fit by a power law in energy 
with photon index 2.40 + 0.11 Jfar ± 0.20 and a flux above 
200GeVof(5.7±0.7 Jffl ,±1.2 JVJ ) x 10~ 12 crrrV 1 . The fit has 
a^ 2 /d.o.f. of 3.5/5. Using the harder image size cut of 200 pe 
used elsewhere in this paper yields consistent results (photon 
index2.29±0.14 J/flf andafluxof (5.5±0.8 JM/ ) x 10" 12 crrrV 1 
above 200 GeV). This flux represents only 2% of the flux from 
the Crab Nebula above the 200 GeV spectral analysis thresh- 
old. Although G 0.9+0. 1 is one of the weakest sources ever de- 
tected at TeV energies, the detection presented here is statisti- 
cally highly significant (13 cr), demonstrating the sensitivity of 
the H.E.S.S. instrument. 




Fig. 2. 90 cm radio flux map of GO.9+0.1 from LaRosa 
et al. (2000 1 (colour scale), overlayed with contours (at 40% 
and 80% peak brightness) of the smoothed and acceptance cor- 
rected count map of gamma-ray candidates (solid lines). The 
H.E.S.S. data are smoothed with a Gaussian of rms 0.03° to re- 
duce statistical fluctuations. The simulated point-spread func- 
tion of the instrument is shown with dotted lines (also 40% 
and 80%). The best fit position of the VHE gamma-ray ex- 
cess is shown with error bars showing combined statistical and 
systematic errors. The innermost circle (dashed) illustrates the 
95% confidence limit on the rms size of the emission region. 

3. Discussion 

X-ray emission from the shell of GO.9+0.1 is only marginally 
detected above the local GC diffuse emission with a flux of 
2 x 10~ 12 erg cirT 2 s in the range 2-10 keV, and is consistent 
with both a thermal or non-thermal origin (Porquet et al. 2003 1. 
The non-thermal flux from the PWN is 5.8 x 10~ 12 erg cm -2 s _1 
in the same energy range, corresponding to a luminosity of ~5 x 
10 34 erg s assuming a distance of 8.5 kpc. The lack of strong 
non-thermal X-ray emission from the shell and the point-like 
nature of the H.E.S.S. detection argue against an origin of the 
VHE emission in the shell. Thus the PWN appears to be a more 
compelling source for the VHE gamma-rays. 

The total power radiated by G 0.9+0. 1 in VHE gamma-rays 
is ~ 2xl0 34 erg s (at 8.5 kpc), compared with ~4xl0 34 erg s _1 
for the Crab Nebula (at 2 kpc) in the energy band 0.2-10 TeV. 
Until now the Crab Nebula represented the only firm identi- 
fication of VHE emission from a nebula powered by a young 
neutron star. 

A simple one-zone inverse Compton model for the VHE 
gamma-ray emission (Khelifi 2002) can be used to explain 
the observed data. A parent population of accelerated elec- 
trons with a broken power-law spectrum is assumed. The mag- 



4 



F. Aharonian et al.: Very high energy gamma rays from the composite SNR G 0.9+0. 1 




10" 1 1 10 

Energy (TeV) 



Fig. 3. Reconstructed VHE gamma-ray spectrum of G 0.9+0. 1 . 
Empty circles show the spectrum derived using the rather tight 
cuts used for figures^andE] Solid symbols show the spectrum 
derived using an image size cut of 80 photoelectrons (optimised 
for maximal significance for 10% Crab strength sources) rather 
than 200 photoelectrons. The line shows a power-law fit to the 
80 pe data set. 

netic field strength B is assumed to be uniform within the 
PWN and is a free parameter. The maximum energy of the 
electron spectrum is kept fixed at 500 TeV as it can not be 
derived due to the absence of measurements in the hard X- 
ray/soft y-ray band. The synchrotron and inverse Compton 
(IC) emission of these electrons are calculated according to 
Blumenthal & Gould dl970t . The Klein-Nishina effect on the 
IC cross section is taken into account. 

For sources in the GC region the presence of radi- 
ation densities considerably higher than typical Galactic 
Plane values leads to substantially enhanced IC emission 
(de Jager et al. 1995 1. The seed photons for IC scattering are the 
cosmic microwave background radiation (CMBR), the galac- 
tic diffuse emission around 100 pm from dust and around 1 pm 
from starlight. The energy spectra of the dust and starlight com- 
ponents are taken from a recent estimation of the interstellar 
radiation field (ISRF) used for the GALPROP code (Strong 
et al. 12000ft . To account for possible local variations of the 
starlight component the energy density of this component is 
kept as a free parameter. The energy density of the dust ISRF is 
kept fixed at 0.23 eV cm 3 as used in GALPROP. The starlight 
ISRF energy density, the electron spectrum parameters and B 
are adjusted such that the computed synchrotron and IC spectra 
match the observations. 

Figure@]shows the broad-band spectral energy distribution 
(SED) of G 0.9+0.1, together with the best model fit. The fit- 
ted mean magnetic field strength is ^6 pG. The model electron 
spectrum has a low energy index of 0.6 and a high energy index 
of 2.9 with a break at 25 GeV. The best fit energy density of the 
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Fig. 4. The spectral energy distribution of the PWN in 
G 0.9+0.1 from radio to VHE gamma-ray. Radio data (trian- 
gles) are taken from (Helfand & Becker 1987J, X-ray data 
(shaded box) from Porquet et al. 2003 The circles show the 
H.E.S.S. data from this work. The solid curve shows a fit of a 
one-zone inverse Compton model to all data. Contributions of 
the CMBR, IR from dust and starlight photon fields to the IC 
emission are shown. 

starlight ISRF is about 5.7eVcirT 3 , that is 50% smaller than 
the value used in GALPROP. This derived value of B is very 
close to the equipartition magnetic field (^5 pG), calculated 
from the whole fitted electron spectrum and given a source of 
radius 1' at 8.5 kpc. Given the observed variation of the X- 
ray spectrum within the PWN, this one zone model is clearly 
somewhat simplistic. However, such a model, describing the 
average properties of the PWN, can describe the available data 
with reasonable physical parameters. Future measurements in 
the GeV regime with GLAST (sensitivity 2 x 10 -12 erg s^'cm" 2 
at 1 GeV) are needed to resolve the ambiguity between the 
magnetic energy density and that of target radiation fields in 
this region. Alternative explanations based on a hadronic ori- 
gin of the VHE emission (for example via the decay of neutral 
pions produced in p-p collisions) are not excluded. 

The detection of VHE gamma-ray emission from G 0.9+0. 1 
provides the first direct evidence for the acceleration of very 
energetic particles in this object. The position and point-like 
nature of the gamma-ray emission, combined with the broad 
band SED of G 0.9+0.1, strongly suggest an origin of the VHE 
emission in the compact central source (previously identified 
as a pulsar wind nebula). 

This detection represents the first step towards the study 
of an emerging population of VHE gamma-ray emitting PWN. 
The detection of new PWN with accurate TeV spectra and with 
morphological measurements will provide information key to 
understanding particle acceleration in the vicinity of pulsars. 
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